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The recently-discovered "ladder" compound LaCu02.5 has 
been found to admit hole doping without altering its struc- 
ture of coupled copper oxide ladders. While susceptibility 
measurements on the parent compound suggest a spin gap 
and a spin-liquid state, NMR results indicate magnetic order 
at low temperatures. These seemingly contradictory results 
may be reconciled if in fact the magnetic state is near the 
crossover from spin liquid to antiferromagnet, and we inves- 
tigate this possibility. From a tight-binding fit to the valence 
LDA bandstructure, we deduce that the strength of the inter- 
ladder hopping term is approximately half that of intraladder 
hopping, showing that the material is three-dimensional in 
character. A mean-field treatment of the insulating magnetic 
state gives a spin-liquid phase whose spin gap decreases with 
increasing interladder coupling, vanishing (signalling a tran- 
sition to the ordered phase) at a value somewhat below that 
obtained for LaCuCh.s. The introduction of an on-site repul- 
sion term, U, to the band scheme causes a transition to an 
antiferromagnetic insulator for rather small but finite values 
of U, reflecting the predominance of (one-dimensional) ladder 
behavior, and an absence of any special nesting features. 

PACS numbers: 75.10.Jm, 75.30. Kz, 75.40. Cx, 75.50.Ee 



I. INTRODUCTION 

One of the interesting and challenging sub-fields of 
low-dimensional quantum magnetism which has emerged 
from the wealth of activity directed at improving the un- 
derstanding of high-temperature superconductors is that 
of ladder systems jy . These consist of n parallel, inter- 
acting chains of S = 5 ions, which can be considered as 
a spin ladder with n legs, and rungs of n-l bonds. The 
ladders have only weak mutual interactions. A combina- 
tion of experimental and theoretical efforts has in the last 
few years produced significant advances in the realization 
and understanding of the properties of spin ladders, sonic 
of which are not at all intuitive. 

Ladder cuprates emerged first with the discovery 
by Hiroi et al. ||] that in the series of materials 
Sr n Cu n +i02n+i it is possible to create two-dimensional, 
stoichiometric copper oxide planes of composition Cu n +i 
02n+i by removing from the uniform CuC>2 plane paral- 
lel, equally-spaced lines of oxygen atoms. It was pointed 
out by Rice et al. || that because these shear defects give 
rise to only weak, ferromagnetic interactions between 
neighboring copper spins, the remaining strips of Cu02 
plane will appear as isolated rt+l-leg ladders of antifer- 



romagnetically coupled spins. These authors proposed 
that the systems should then illustrate the contrasting, 
and now well-established properties of even- and odd-leg 
ladders, that the former show a gap to spin excitations 
(spin gap) with consequent exponential spatial decay of 
correlations, while the latter are gapless with power-law 
decays. Subsequent susceptibility Q and nuclear mag- 
netic resonance (NMR) || experiments have amply borne 
out this conjecture. 

The theoretical understanding of ladders has been in 
large part based on numerical techniques, which are par- 
ticularly well suited to systems of such restricted dimen- 
sionality, and this is reviewed in Ref. pj. The first in- 
dications that the 2-leg ladder should exhibit a spin gap 
came from numerical work HQ] , and a variety of meth- 
ods has since been applied to investigate the spin gap and 
spin correlations, and the properties of multi-leg ladders. 
|p|,ppC|] These studies not only confirm the picture of the 
spin system emerging from analytic approaches valid 
in certain limits, but provide the most accurate informa- 
tion available on the properties of this class of strongly- 
correlated systems. 

The second proposal concerning the properties of lad- 
der compounds made in Ref. |3|, that the doped lad- 
der should become superconducting with a e?-wave order 
parameter, has proved harder to test for materials rea- 
sons. However, the compound LaCu02.5, recently syn- 
thesized by Hiroi and Takano, jl2j has been found to 
admit hole doping without altering its structure of cou- 
pled, 2-leg copper oxide ladders, and so constitutes the 
first case in which one may seek doping-dependent behav- 
ior analogous to that of the high-temperature supercon- 
ductors. The high-pressure phase of Lai_ x Sr :c Cu02.5 is 
derived from a cubic, three-dimensional perovskite struc- 
ture. The absence of oxygen atoms along lines leaves 2- 
leg ladders which relax from a relative angle of 90° in 
the primitive structure to 62° in the depleted one. The 
material is then orthorhombic (space group Pbam), with 
four copper atoms per unit cell, and may be considered 
as a set of ladders oriented along the z-axis, and peri- 
odically arrayed in the x- and y-directions, as shown in 
Fig. 1(a). This structure is not altered with Sr doping 
to x = 0.2. The interladder coupling in the (x,y) plane 
arises because each oxygen atom at the outside edges of 
the 2-leg ladders, which is part of the planar CUO4 unit, 
is also effectively apical to a copper atom in a neighboring 
ladder, contributing to a finite transfer integral. 

Susceptibility measurements [fL2| on the parent com- 
pound up to temperatures of 500K interpreted by a for- 
mula proposed by Troyer et al., [fL3| suggest the pres- 
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ence of a spin gap in the excitation spectrum, and there- 
fore a spin liquid state. In contrast, NMR studies of the 
same samples indicate that the system orders antifer- 
romagnetically at low temperatures, below an apparent 
Tjv — 117K. Since the spin susceptibility should evolve 
continuously as one passes through the quantum critical 
point separating the spin-liquid state and the state with 
long-range antiferromagnetic order, it is possible that 
even in the latter, close to the quantum critical point, 
the spin susceptibility will decrease substantially as the 
temperature is lowered. Such behavior may be difficult 
to distinguish from that of a true spin liquid with a finite 
gap. Therefore it is plausible that the seemingly con- 
tradictory experimental results may be reconciled if the 
magnetic state is close to the crossover from an antifcr- 
romagnet to the spin liquid. Here we seek evidence, by 
examining the effects of interladder interactions, that the 
system is indeed close to this quantum critical point. 

The outline of this paper is as follows. In section II 
we present a tight-binding fit to the LDA bandstruc- 
ture to extract the interladder hopping matrix elements, 
and use these to obtain the superexchange interactions. 
We consider in section III the nature of the spin-liquid 
ground state for ladders of spins coupled in both or- 
thogonal directions, using a mean-field approach to es- 
timate the location of the quantum critical point. In 
section IV we introduce a double-occupancy term to the 
tight-binding bands, to investigate the degree of three- 
dimensional character in the electronic structure as a re- 
sult of interladder interactions. Section V contains our 
conclusions and a brief discussion. 



II. TIGHT-BINDING FIT TO THE LDA 
BANDSTRUCTURE 

The bandstructure of LaCuC>2.5 has been computed 
by Mattheiss, jl5| using the Local Density Approxima- 
tion (LDA) method. Here we examine the dispersion of 
only the highest occupied valence bands, and use a tight- 
binding model based on the single, copper-based orbital 
in each planar CUO4 unit within the ladder which lies 
closest to the chemical potential: in the Cu n +i02n+i 
strip this is the antibonding, Cu-centered 3d x 2 _ y i or- 
bital. Restricting the set of hopping matrix elements to 
those between nearby copper atoms, we obtain an effec- 
tive one-band model, albeit with four mixed levels (one 
from each Cu atom in the unit cell). These generate a 
complex of four bands, which lie close to or cut the Fermi 
energy. The complex is half-filled in the undoped system, 
and well separated from other bands on the small energy 
scales of most physical interest, so can be taken to de- 
termine the low-energy behavior of the model. With this 
interpretation, one may deduce the ratio i'/i, of inter- 
and intraladder atomic orbital overlap, and thus estimate 
the ratio J'/ J of the magnetic interactions, by using the 
superexchange result [jl6| J ~ At 2 /U. 

The tight-binding Hamiltonian is 



H 



(i) 



in which i and j each denote a pair (n, m), where n labels 
the unit cell, and m = 1, ... 4 the different atoms within 
each cell. The hopping matrix elements Uj may be taken 
to be short-ranged, and the maximal set which we need 
to achieve a reasonable fit is illustrated in Fig. 1(b). 
Starting from tight-binding parameter fits in the two- 
dimensional CuC>2 plane, |{l7|,[l8) we choose the nearest- 
neighbor intraladder parameters, t r for hopping along a 
rung and t z for hopping along a leg, to have the same 
value t r = t z = t ~ 0.4eV. We use as the next-nearest- 
neighbor parameters t' r — —\t r and t' z — \t z for hop- 
ping across a square plaquette in the ladder, and along 
two leg bonds, respectively. Finally, we introduce a pa- 
rameter t s for hopping from a copper atom on one lad- 
der to its neighbor on the adjacent ladder, and also the 
next-neighbor analog t' s for transfer between atoms on 
adjacent ladders with a relative displacement of one leg 
bond. Because of the low symmetry of the LaCuC>2.5 
structure, the Cu-(apical) O-Cu interladder bond is far 
from straight, with the CUO5 pyramids quite irregular, 
and this distortion may allow a substantial value of t s . 
The role of each of the terms in the fitting scheme will 
be illustrated below. The next-neighbor parameters t' r a z 
are expected to be significantly smaller than their direct 
counterparts, due to the short-range nature of the over- 
laps, and the important parameter to fit will be the ratio 
of t s to t r . 

The Hamiltonian ([!]) may be expressed in matrix form 
in reciprocal space as 



H = 4 H k c ko 



(2) 



where cl = (c k t., c?L, , c^t.) is the vector of creation 
operators 



mt 
C kJ 



(3) 



for the Bloch states formed by separate linear combina- 
tion of the d x 2_ y 2 orbitals of each copper atom in the 
unit cell, and k is a vector in the orthorhombic Brillouin 
zone. Setting the lattice constants a, b and c to unity, 
the Hamiltonian matrix is 



( t z (k z ) t s e lk y>- 

ts Q-ik v s y t z (k z 



t r e~ ik - r 



\ t r e 





zk.r 



t r e 





t r e lkr 

_t z (k z ) t s e* k y s y 

tse -ik ySy tz ( k ^ j 



(4) 



where t z (k z ) = 2t z cos k z + 2t' z cos 2k z , t s = 2t s cos^k x , in 
which the factor cos^k x arises because s x — 0.5, t v = i„ + 
2^cosfc z , and the bond vectors r(r) and s(s) are shown 
in Fig. 1(c). The eigenvalue problem gives an equation 
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quadratic in the squares of the mode frequencies, whose 
solutions are the dispersion relations of the four energy 
bands 

e k = ± [tl + At 2 . cos 2 \k x ± At r t s cos \k x cos ±k y ] ^ 2 

— 2t z cos k z — 2t' z cos 2k z . (5) 

This concise, closed form emerges because the exponen- 
tial factors e M 2ik v r v+ 2ik v s v) ma y be collected as 1— cos k y , 
as r y + Sy = 0.5. The general dispersion simplifies further 
for particular values of k x and k y in the Brillouin zone: in 
particular, on the zone faces (k x ,k y = tt), the last term 
in the square root vanishes and the bands are doubly de- 
generate, as required by the group-theoretical analysis of 
structures with non-symmorphic space groups. 

The dispersions of the 4 bands are shown in Fig. 2(a) 
for a series of high-symmetry lines in the orthorhombic 
Brillouin zone, and for the ratio t s /t r — 0.4. The la- 
belling of points is shown in Fig. 2(c), and their order 
is chosen to match the results of Matthciss. Also 
shown in the figure are the energy-spacings dictated by 
the choice of the tight-binding parameters, as these vary 
between their maximum and minimum values along the 
chosen directions. 

For an isolated ladder along z, i.e. with no interladder 
interactions (t s = 0), the bands would be completely flat 
around TXSYr and ZURTZ, with only a cosine disper- 
sion (mildly perturbed by the t' z term) along TZ. In this 
situation there would be two doubly-degenerate bands 
around the zone center and zone face, corresponding to 
the bonding and antibonding bands of each ladder. These 
features are reflected in the two branches observed in Fig. 
2(a). From Eq. (^), the separation of the degenerate 
bands along the zone face XSY and edge URT is gov- 
erned by the combinations t r ± 2t' r , and the splitting of 
the degeneracy along TX, TY and ZU, ZT by the combi- 
nations t s ± 2t' s . Considering first the intraladder param- 
eters t z , t' z , t r and t' r chosen from the two-dimensional 
CuC>2 plane, we find good agreement of the tight-binding 
result with that from LDA. Particularly notable is that 
the negative sign of t' r , and the relative magnitude jrt r , 
are required to reproduce the bands separations along 
XSY and URT. There is no evidence that a value of t r dif- 
ferent from t z would improve the fit. The next-neighbor 
hopping parameter along the ladder, t' z , appears only as 
an asymmetry of the cosine dispersion along TZ, and the 
value chosen is in qualitative accord with the LDA result. 
That Cu02 plane parameters remain appropriate for the 
ladder confirms the predominantly local picture of the 
interactions between copper sites. Turning to the band 
splitting at the T and Z points, we find that the relatively 
large value t s = \t r (Fig. 2(b)) gives the best qualitative 
reproduction of the bands crossing the chemical poten- 
tial in the full LDA calculation by Mattheiss, |15| but 
that the value t s = \t r appears closer to the results of a 
"I2-parameter fit" illustrated in the same reference. The 
difference of the T and Z splittings is given rather well 
by the starting choice of t' s — ^t s . 



It is clear from Fig. 2 that the primary feature of the 
dispersion remains that in the k z direction, i.e. along the 
ladder. The Fermi surfaces for each band are determined 
almost exclusively by this part of the dispersion, in that 
they appear as sheets which are almost flat, perturbed 
little by the t r and t s terms, and have k z as normal. At 
half-filling, the lowest band is an exception to this situa- 
tion, because part of it is also filled close to the Z point. 
On doping with holes, this region is rapidly emptied (be- 
low 5%) so that all four Fermi surfaces are sheets parallel 
to the (k x , k y ) plane. Only when the doping level reaches 
20% does the chemical potential drop below the highest 
band in the TXSYr plane, causing a pocket to open at 
the r point for this band. 

In summary, we find that a simple model of a sin- 
gle orbital per copper atom provides a good fit to the 
bandstructure. While there is scope for some variation 
in the choice of intraladder parameters, with this level of 
agreement between the tight-binding results and those of 
LDA it is not worthwhile to optimize further. We choose 
to work with the above values of the intraladder hop- 
ping matrix elements, and with the interladder overlap 
ts = \tr, bearing in mind that this latter will be close to 
the upper limit of the narrow range of probable values. 
Estimating the superexchange interaction by J' oc At' 2 /U 
leads us to conclude that the interaction between spins 
on neighboring ladders will have a magnitude J' ~ 0.25 J, 
where J is the intraladder magnetic coupling of both rung 
and leg spins. Because J' is an appreciable fraction of J, 
it is clear that the spin interactions in LaCu02.5 will have 
significant three-dimensional character. 

III. MEAN-FIELD ANALYSIS OF THE SPIN 
GROUND STATE 

A mean-field analysis of the spin state for ladder sys- 
tems was introduced by Gopalan et al. [p~9|| , and in this 
section we follow closely the treatment of these authors. 
They employed a bond operator representation of S = h 
quantum spins, used initially ppj ] to investigate dimer- 
ized spin phases in two-dimensional systems, exploiting 
the fact that the topology of the ladder favors dimeriza- 
tion, particularly when the spin interaction, J, on a rung 
exceeds that on a leg, A J. The authors proceeded in the 
mean-field approximation to consider first the properties 
of an isolated ladder as a function of the interaction ratio 
A, then of periodic arrays of ladders in two dimensions, 
and finally of an array of frustrated double ladders of the 
type found in the SrCu 2 3 system. [|J This approach can 
be considered to be exact in the limit where A — > and 
the spins form dimers on every rung, while its accuracy 
will diminish on extrapolating through finite A towards 
the desired isotropic point, A = I. 

We begin by representing the system of spin ladders, 
each coupled by two identical bonds per spin to sepa- 
rate, neighboring ladders, in a geometrically simpler but 
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topologically identical form, as shown in Fig. 3. Note 
that the spin configuration for local, antiferromagnetic 
interactions is unfrustrated, so a transition to an ordered 
antiferromagnet may be expected with increasing inter- 
ladder coupling. The magnetic interactions are taken to 
be J for spins on the same rung, A J for spins separated 
by a leg bond, and A' J for neighboring spins on different 
ladders, while no other couplings are considered. The 
Hamiltonian for the spins is 

H = J S Si j.S r j + XJ S S m j.S m j_)_2 

+A'jWs rj -.S i .,i. + S r j.S ...1. i.V (6) 



2 2 • 



where j is a rung bond index and the labels I and r 
denote spins on the left and right sides of the ladder. 
Following Ref. Jl9| , transformation to the bond-operator 
representation yields 



where 



H — H + Hi + H 2 + Hho, 

H ° = J E(-! s ] s i+i*L^) 
-E^ ( s i s J ~ 1 



(7) 



(8) 



(9) 



and 



i i \ -a J • 2 ^ '2 ! 

In these equations, st is the creation operator for a spin 

singlet on bond j, the operators ft Q create the three pos- 
sible triplet states on the same bond, and the Lagrange 
multiplier, fij, introduced to ensure the constraint 



S l S 3 n ~ / './■'■' 



(11) 



on each bond, which restricts the physical spin states to 
singlets or triplets, appears as an effective chemical po- 
tential. The part Hj w in (Q) contains terms with three 
and four tj >a operators, and will be neglected in our ap- 
proximation; in Ref. [19J it was shown in addition that 
the effects of such higher-order terms are small. 



Because the singlet on each bond has the lowest energy, 
we assume that the system condenses into this state, lead- 
ing to a finite expectation value of the bosonic Sj oper- 
ator, (sj) = s. This is the average expectation value, 
or mean- field value of the operators Sj, and the site- 
dependent chemical potential fij is also replaced by a 
global average value \i. Working with the physical unit 
cell (Fig. 3), which contains two rungs, we may trans- 
form the operators tj iCt in Hq, Hi and Hi to those for 
two types of triplets, t\ and t 2 . The Hamiltonian in 
this approximation is [£lj 

H m (li, s) = N {-ps 2 - M* 2 + M) + E | E [ A k*kl*L 



ka Kv=l,2 



in which 



Ak = | J - fi + Js 2 Acos k z , 



Ak = \ Js 2 Xcosk z , 



A k = 2A^ = — J's 2 cos \k x cos \k y , 



[K^Av (12) 

[1-2]}, 



(13) 



(14) 



(15) 



and N denotes the total number of ladder rungs. The 
part of H m (|l2|) dependent on the triplet operators is di- 
agonalized by the non-unitary, bosonic Bogoliubov trans- 
formation 



lL = cosh 9£ (ti a ±tl a )+ sinh o£ ( ± t x X a - t\ a 



(16) 

whose coefficients are given by 
cosh2^ = ^4^ sinh2^ = ^4^), (17) 



where in turn 



(A k ±A^) 2 -4(A k ±A^r 



(18) 



yield the dispersion relations of the two magnon 
branches. The Hamiltonian now takes the form 

H m ( f x,s)=N (-§ Js 2 - (is 2 + (i) - |iV {\J — n) 

+ EE^(^7L + |), (19) 

ka f=± 

which contains the mean-field part and the zero-point 
quantum corrections from the triplet magnon excitations. 
The mean-field equations to be solved self-consistently for 
fi and s are given by 



4 



and 
,dH« 



5 x 

k v=± 



(A k + »K) 



n m K) (20) 



k i/= ± 



Ak - 2A t 
2< 



a k n »n( w k)> 



where 



= A cos fc z ± A' 



COS 2" ^2; 



(21) 



(22) 



contain the dispersive parts of ui^. n m {ui^) denotes the 
magnon thermal occupation function, and will be dis- 
cussed in more detail in a future publication. The factor 
of 3 preceding the k-summations in both equations is 
the result of the sum over a for the three triplet magnon 
states. p3] This factor was omitted in Ref. |^9| , and we 
comment below on the effect of the correction on the re- 
sults presented there. 

The mean-field equations are solved at zero temper- 
ature, where the thermal factor becomes unity, and by 
taking the continuum limit in which the k-sum becomes 
an integral over three-dimensional reciprocal space. As 
in Ref. we reduce the two equations to a single one 
for the variable 



d = 



2Js 



t:2 



which has the form 



d = 5-f 



E 

v=± 



d 3 k 



(2tt) 3 ^1 + dal 



(23) 



(24) 



As a characteristic parameter of the spin-liquid ground 
state, we will be most interested in the value of the spin 
gap, the minimum excitation energy of the lower branch 
of the triplet magnon excitations, which is given by 



A = (i J - fl) y/l-d(\ + \>), 



(25) 



with d determined by the mean-field equation (|2J) . From 
Eq. ( |l8| ) we see that the excitation spectrum has a min- 
imum at the wavevector kjf = (0,0, tt) in the reciprocal 
lattice of the bipartite structure shown in Fig. 3. The 
value of A' where A is driven to zero will give the transi- 
tion from the spin liquid state, where the spin orientation 
fluctuates with a time-averaged value of zero and with 
short-range correlations primarily at kj/, to a magneti- 
cally ordered state characterized by kjvf . This wavevector 
corresponds to uniform polarization of the spin singlets 
on ladder rungs in the (x, y) plane, with spins oppositely 
directed between neighboring planes in the z-direction, 
i.e. a simple antiferromagnetically aligned spin pattern. 

In the limit of no interladder coupling we obtain the 
spin gap A = 0.501J of the isolated, isotropic (A = 1) 2- 
leg ladder. This is only a mean-field result, but is in very 



good agreement with the result Ao = 0.504 J of numerical 
studies || by the Density Matrix Renormalization Group 
technique. In fact this agreement is largely serendipitous, 
and deteriorates on taking into account the higher-order 
terms p2 ] ; in the mean-field approximation, the spin gap 
of the isolated ladder diverges logarithmically in the limit 
of large A, |l9| and the effects of this increase are already 
manifest at A = 1, causing the mean- field result, which 
initially underestimates the spin gap, to recover towards 
the exact value determined numerically. In this treat- 
ment we have not had to invoke a self-energy correction 
term: in Ref. Jl9[ , the authors investigated the curious 
qualitative behavior of the solution to their (erroneous) 
mean-field equation by expanding in small A about the 
limit of strong rung coupling where the dimer treatment 
is accurate. In the corrected mean-field theory, one ob- 
tains 

A = J(l-A+ ±A 2 + fA 3 + 0(A 4 )), (26) 

which corresponds reasonably well to the result 

A = J(l-A + ±A 2 + ±A 3 + 0(A 4 )) (27) 

of a detailed strong-coupling analysis including excitation 
modes. 
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Previously, the coefficient of the quadratic 
term had been found to be negative, and so a self-energy 
term [3\ 2 was introduced to correct for short-range in- 
teraction effects which appeared to have been missing at 
the mean-field level; the chosen value (3 — 0.7 brought 
the results into good argeement with previous numerical 
ones, and with the above approximate treatments. 

The spin gap obtained from the solution of (53) for 
the three-dimensionally coupled ladder system is shown 
in Fig. 4(a) as a function of the interladder coupling A', 
for fixed A = 1. We see immediately that the spin gap 
decreases monotonically, with the transition point at A' = 
0.121. Comparison with the result of section II indicates 
that the LaCuC>2.5 system should lie within the ordered 
antiferromagnetic regime, but that it is indeed located 
in the vicinity of the quantum critical point marking the 
phase transition from spin liquid to magnetic order. 

It is also instructive to compare the appearance of the 
spin gap with that in a two-dimensional, unfrustrated, 
periodic array of coupled ladders. |nj In this case there 
is only one magnon branch, and the dispersive factor in 
the excitation spectrum is 



«k = A cos fez 



5 A' cos k x 



(28) 



where k x is a wavevector parallel to a ladder rung. Again 
kju = (0, w) corresponds to a simple antiferromagnetic 
spin alignment in the ordered phase, and the factor of 
i appears because there is only one bond per spin be- 
tween ladders in such a system. The spin gap for the 
two-dimensional array is shown as a function of A' in 
Fig. 4(b), where again A = 1 and the gap at A' = is 
that of the isolated, isotropic ladder. The critical value of 
the interladder coupling, \' c — 0.43, is seen to be signif- 
icantly greater than twice that in the three-dimensional 
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case above, as might be expected for the simple rea- 
son that there are half as many interladder interactions, 
which may be taken as an indication that the spin liquid 
state is more robust in lower dimensions. The almost 
linear decrease, in contrast to the downward curvature of 
the function in Fig. 4(a), illustrates a further effect of di- 
mensionality, and agrees well with the qualitative result 
of Ref. || . 

We may conclude that to within the accuracy of the 
mean-field approach for an isotropic ladder system, the 
above analysis provides good evidence that LaCu02.5 nes 
on the magnetically ordered side of the quantum critical 
regime of the transition between spin liquid and antifer- 
romagnet. 



atoms in the two (x, y) planes of the doubled unit cell, 
and 



H 



k = 



/ M+ S(k z ) R (k z ) \ 

S*(-k z ) M_ R(k z ) 

R*{-k z ) M+ S(k z ) 

V TL(-k z ) S*{-k z ) M_ J 



in which the 2x2 matrices 



M± = 



It' cosk z ± W5n t z e l 2 



2t' z cosk z =F hUSn 



(30) 



(31) 



IV. HARTREE-FOCK APPROXIMATION TO 
THE HUBBARD MODEL 

We have seen in the previous sections that the interlad- 
der interactions in LaCuC>2.5 are quite strong, and possi- 
bly sufficiently strong to change the magnetic structure 
from a spin liquid to a three-dimensionally ordered anti- 
ferromagnet. To examine further the character of the 
electronic structure, we introduce an on-site Coulomb 
interaction, U, to the tight-binding bandstructure of 
section II, and perform a Hartree-Fock calculation of 
an ordered antiferromagnetic state. In a strictly one- 
dimensional system, the critical value of U necessary to 
stabilize an insulating antiferromagnet vanishes at half- 
filling, whereas in a general, three-dimensional electronic 
structure it is of the order of the bandwidth. 

Introduction of the on-site Coulomb interaction leads 
to a Hubbard Hamiltonian 



HM 



ija 



(29) 



Examining the stability of an antiferromagnetic state 
with wavevector \<lm , we note first that there will be no 
increase of the unit cell in the (x, y) plane, since there 
are already two atoms for each spin direction, but that 
it doubles along z when the spins order along the lad- 
der legs. In the notation of A and B sublattices for the 
bipartite system, sites m = 1,3 of the original unit cell 
(Fig. 1(c)) in every second plane, and m = 2,4 in the 
alternating planes, will belong to the A sublattice, while 
the remaining sites will belong to B. Introducing the pa- 
rameter Sua — Hf — 7i j as the difference between average 
site occupation by particles of each spin orientation on 
a site of the A sublattice, we require Sua = —Sub = Sn 
for all sites. We proceed by solving the problem in the 
Hartree-Fock approximation, with the value of U where 
Sn becomes finite marking the antiferromagnetic transi- 
tion. 

Following the treatment of section II, the Hamilto- 
nian may be written as in Eq. 



with now 



„M l&t 



R(fe.) = 



Ue 1 ^ t'e zk - r+i 2 fe * 



tie 



ik.r 



(32) 



R(k z ) (defined identically using r) and 
/ t s cos ±k x e lkySy 



S(k z 



t' cos±k x e lkySy+i 2 k * 



t s cos ^kx^ v v 



l 2 kz 



t s cos ^k x e lkySy 



(33) 



are the generalizations of the previous expressions to the 
new unit cell, and the new c-axis dimension is set to unity. 
This matrix cannot be block-diagonalized, but the struc- 
ture of the solution for the eigenmodes and eigenvectors 
is evident from section II. Schematically, if the (k x ,k y ) 
dispersion contained in the square root in Eq. (||) is de- 
noted by e£ , the eight band dispersions will have the form 

El = ± [el + U 2 Z cos 2 \k z ± \U 2 8n 2 } 1/2 - 2t' z cos 2k z , 

(34) 

where i labels the bands. The U Sn term splits the former 
four bands into two sets, which as U becomes large will 
not overlap, ensuring that the half-filled system becomes 
insulating. The corresponding eigenvectors are the states 



(35) 



where the index j runs over the eight atomic sites, and 
is the Bloch state created by the operator c£. The 
Hartree-Fock equations for the system, which will deter- 
mine 5n and the chemical potential fj, at fixed U when 
solved self-consistently, are the equations for the total 
and site occupancies 



1 



\ ^ 



d 3 k 



1 



I — 1 



(36) 



where the superscripts a and b denote an( j 
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8 r rl^h \CU 2 
Kl-flP + faJ.g/^-lgit-. (37) 

These equations have been generalized to arbitrary band 
filling, which is expressed in terms of the deviation S 
from half- filling, which in turn is normalized to be 1. 
The chemical potential is contained in = — M) an d 
the second expression is valid for the coefficients of any 
chosen Bloch function j. 

The operation of diagonalizing the Hamiltonian ma- 
trix (^) can be performed numerically at sufficient speed 
that it is still possible to solve the Hartree-Fock equa- 
tions, which involve three-dimensional k-integration, on 
a workstation within a reasonable amount of time. In 
Fig. 5(a) are shown the energy bands for the incipient an- 
tiferromagnetic system for half-filling (<5 = 0) and U = t. 
The splitting of the bands into upper and lower branches 
as a result of U is clearly evident, as is the characteris- 
tic TXSYr and ZURTZ structure of the purely kinetic 
Hamiltonian in both sets of bands, at positive and nega- 
tive energies, as a result of the folding back of bands due 
to the change in meaning of the coordinate k z , which now 
spans a Brillouin zone half the former size. At this value 
of U, the half-filled system would appear to be close to 
the transition from metallic to insulating behavior, which 
one expects near the point where there is no longer any 
overlap between upper and lower band energies in any 
region of reciprocal space. At smaller values of U, the 
two sets of bands are characteristic of the doubled Bril- 
louin zone, and have a semi-metallic overlap, while as U 
is raised to large values, the energy gap increases, and 
the bands become progressively more flat. 

In other cuprate compounds, a Hubbard model has 
been found to give an accurate description of the 
low-energy behavior with a ratio U/W of order unity, 
where W is the bandwidth. Here the total bandwidth of 
the uppermost valence bands is approximately the same 
as in the cuprates with CuC>2 planes, W ~ 3eV. Since the 
local environment of the Cu 2+ ions is similar, the on-site 
Coulomb repulsion should also be the same, U ~ 4eV. 
These parameter values place LaCu02.5 we U within the 
Mott insulating region. We note that the actual magnetic 
structure cannot be determined in the Hartree-Fock Ap- 
proximation, as the quantum corrections which act to 
stabilize the spin liquid phase are not included. 

The form of the Hartree-Fock solutions for smaller val- 
ues of U are sensitive to the effective dimensionality of 
the electronic structure. In general, values of U ~ W are 
required to obtain a Mott insulator at half-filling, but 
the perfect-nesting property of a one-dimensional band 
gives an insulating state for arbitrarily small values of 
U. In Fig. 5(b), the average antiferromagnetic order 
(parameterized by Sn) in the half-filled system is shown 
as a function of the on-site repulsion U, which is mea- 
sured in terms of the bandwidth W = 3eV. In LaCuC>2.5 
we find that the critical value U c for the antiferromag- 
netic transition is rather small, U c ~ 0.2IF, which in- 
dicates substantial one-dimensional nature. The inter- 



ladder hopping matrix elements, although quite strong, 
are not sufficient to destroy the nesting character com- 
pletely, which would make the electronic structure effec- 
tively three-dimensional. 

This quasi-one dimensional behavior of the bands may 
also be reflected in the sensitivity of the system to ran- 
dom potential fluctuations. It is well known that the 
onset of localization is strongly dependent on dimension- 
ality. Clearly, in ladder compounds there is an inherent 
conflict between the need to change the valence of the 
countcrions in order to induce hole carriers, and the need 
to avoid strong, random potential fluctuations. Thus it 
appears that in Lai_ a; Sr a ;Cu02.5 the random potential 
fluctuations act to cause localization for x < 0.15, in 
spite of the substantial interladder overlap, suggesting 
that a more gentle hole-doping technique will be required 
to retain itinerant character at small doping concentra- 
tions. This could perhaps be achieved in structures where 
the counterions are further from the ladders than in the 
present case of LaCu02.5- 

V. CONCLUSION 

We have investigated the basic electronic and magnetic 
properties of the three-dimensionally coupled 2-leg lad- 
der compound LaCu02.5, a material which is of signifi- 
cant experimental and theoretical interest as it is the first 
ladder compound to be discovered in which the CU2O3 
ladders may be doped with holes. We present a tight- 
binding fit in which the bands are derived from a single 
(d x 2_ y 2) orbital close to the Fermi energy on each copper 
atom. The results of LDA studies are well reproduced by 
reasonable values of the most significant transfer inte- 
grals: those within each ladder are found to be similar 
to the Cu02 planar system, emphasizing the short-range 
nature of the dominant physical processes, while the in- 
terladder hopping term t' ht is found to be quite large. 
As a consequence, the interladder magnetic coupling J' is 
also relatively large, and the compound may be expected 
to exhibit some three-dimensional characteristics. 

The effective spin interactions in this structure are 
those of unfrustratcd antifcrromagnetism, and a mean- 
field treatment of the magnetic state from the basis of 
dimerized singlets on the rungs of decoupled ladders gives 
a spin-liquid phase whose spin gap decreases with increas- 
ing interladder coupling. The spin gap is found to vanish, 
signalling a transition to an ordered phase, at an inter- 
ladder coupling ratio J' j J somewhat smaller than that 
deduced for LaCu02.5, indicating that the system is lo- 
cated in the antifcrromagnetically ordered state, albeit 
not far from the quantum critical point of the ordering 
transition. We may take the very small value of the in- 
trinsic susceptibility measured at low temperatures 
as evidence that proximity to the critical point, and the 
possibility this allows of significant critical fluctuations, 
plays an important role in determining the physics of the 
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system. 

Further insight into the electronic properties is pro- 
vided by the introduction of an on-site repulsion term, 
U, to the band scheme: within the Hartree-Fock approx- 
imation we find that in the half-filled system the transi- 
tion to an antiferromagnctic insulator occurs for values 
of U quite small compared to the bandwidth. This is 
a reflection of the fact that the predominant feature of 
the bands remains the dispersion in the ladder direction, 
and this one-dimensionality, also apparent in the Fermi 
surfaces of the partially-filled bands, makes the system 
inherently susceptible to potential fluctuations. That U c 
is finite illustrates an absence of perfect nesting condi- 
tions in the bandstructure, and the instability occurs at 
the antiferromagnetic wavevector. 

Our results represent the first theoretical consideration 
of the LaCu02.5 system. While by the nature of the ap- 
proximate techniques used they are somewhat inexact, 
we believe that they are important in establishing the 
parameter space for the magnetic state, and in focussing 
the direction of further research. The properties of such 
"nearly critical" magnetic systems have been studied for 
the two-dimensional, planar case by detailed analytical 
PH and numerical |Q techniques. We propose the ap- 
plication of similar methods for the case of ladder systems 
with variable interladder coupling in two additional di- 
mensions, to investigate the nature of the critical point, 
and the surrounding "quantum critical" regime, in the 
same formalism, with particular emphasis on the appear- 
ance of physical quantities accessible in experiment. Such 
studies may prove valuable in establishing the framework 
within which to interpret the data from experiment, and 
hopefully will serve to reconcile the apparent contradic- 
tion in current results. 

The observation that the spin liquid is not the ap- 
propriate description of the magnetic ground state in 
LaCu02.5 is itself important, particularly with regard to 
the current interpretation of the susceptibility data. ]T^ ] 
While it may be possible to deduce the form of the sus- 
ceptibility suitable for the critical regime from a mean- 
field picture, we await the results of numerical studies of 
the same system in order to make a more detailed com- 
parison with the data. A further direction not strongly 
emphasized in this communication is the nature of the 
doped system: the methods used here become less accu- 
rate at finite doping, and so were not studied in great de- 
tail in this regime, but may nonetheless be used to obtain 
additional insight into the evolution of spin properties on 
moving towards the metallic state. 

We are grateful to Y. Kitaoka, M. Sigrist, M. Takano, 
M. Troyer and especially Z. Hiroi and M. Zhitomirsky for 
helpful discussions, and to the Swiss National Fund for 
financial support. 
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FIGURES 



FIG. 1. (a) Crystal stucture of the depleted perovskite 
LaCu02.s- Black and white spheres within the ladder units 
represent Cu and O atoms, respectively, and grey spheres 
represent Sr. The material is viewed along the axis of the 
ladders (z), which can be seen to be rotated about this axis 
to a relative angle of 62°. In any (x,y) plane each cop- 
per atom is bonded to a neighbor in the same ladder by 
a rung bond, and to two copper atoms in neighboring lad- 
ders. The interladder couplings are through one oxygen atom 
which is bonded as part of the square planar coordination in 
the same ladder and is apical to a copper atom in the next 
ladder, and through the single apical oxygen atom, which 
is bonded in the ladder of the other neighbor. Allocation 
of antiferromagnetically arranged spins to each ladder shows 
that the material can exist as an unfrustrated antiferromag- 
net. (b) Schematic representation of LaCu02.5 to show the 
tight-binding parameters between Cu atoms used to fit the 
LDA bandstructure. (c) Appearance of the four inequiva- 
lent Cu atoms (black circles) in the unit cell. White circles 
represent O atoms. The vectors for the two types of bond 
in the (x,y) plane are r(r) = (0.5856a, -(+)0.2114fe) and 
s(s) = (+(-)0.5a, 0.28866), where a ~ y/2a p and b ~ 2y/2a p 
are the lattice constants in the x and y directions, and a p — c 
is the lattice constant of the original, cubic perovskite struc- 
ture, n 



FIG. 5. (a) Bandstructure computed in the Hartree-Fock 
approximation for the parameters of the tight-binding fit of 
section II, with the inclusion of an on-site repulsion param- 
eter U. The bands are shown for U — t — 0.4eV, where the 
half-filled system is close the metal-insulator transition, (b) 
Antiferromagnetic order as a function of U. The degree of 
ordering is parameterized by 8n — n-\ — nj, , the difference be- 
tween the average occupation of each site in the structure by 
particles with spins oriented upwards and downwards. The 
doping is S = 0, i.e. the system is half-filled, and the temper- 
ature is taken to be low. 



FIG. 2. (a) Illustration of the effects of the chosen band- 
structure parameters on the observed dispersion curves, for 
the choices t z = t r = t = 0.4eV, t' z = t' r = — ^t r , t s = |t r 
and t' g =\t s . (b) Tight-binding bandstructure for the param- 
eter set which appears closest to the LDA results of Mattheiss. 
JlB"! Parameters are as in (a), but with t s — \t r - Details of 
the fitting procedure and parameter choices are given in the 
text, (c) Notation for k points in the orthorhombic Brillouin 
zone. 



FIG. 3. Schematic representation of the periodic structure 
in the (x,y) plane of locally antiferromagnetically-correlated 
spins, for calculation of the spinon dispersion. The magnetic 
interaction parameters shown are J on the ladder rungs and 
J' = A'J between spins on neighboring ladders. 



FIG. 4. (a) Spin gap A as a function of the ratio be- 
tween interladder and intraladder magnetic couplings J' /J, 
calculated for a three-dimensional system at T = 0, in the 
mean-field approximation from the starting point of dimerized 
ladder rungs. Here the intraladder rung and leg interactions 
are the same (A = 1), and the value of the spin gap at J' — 
is that of the noninteracting, 2-leg ladder (Ao = 0.501J) in 
the same approximation, (b) Spin gap for a two-dimensional 
array of ladders. 
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